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S^' Abstract 

We explore the prospects of observing leptonic CP violation in a neutrino 
^ ' factory in the context of a scenario with three strongly oscillating neutrinos able 

to account for the solar, the atmospheric and the LSND results. We address also 
the problems related with the fake asymmetries induced by the experimental 
device and by the presence of matter. 
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1 Measuring leptonic CP violation with neutrino 
beams from a muon collider 



Muon storage rings, muon colliders, and their physics potential are being studied inten- 
sively at FNAL and at CERN [|l|. In particular, a muon storage ring at some 20 GeV, 
being the first step in these projects, would serve to produce intense neutrino beams 
of unique quality. This possibility has received much attention p|, |^, |^, ^ recently. 
The straight sections of a muon collider would serve as sources of I/^ and of v^. with 
energy spectra perfectly calculable from muon decay, when positive muons are stored, 
and, similarly, of z/^ and V^. with well-known energy spectra when negative muons are 
stored. 

In this work we explore the possibilities of studying CP violation in the leptonic 
sector, at a neutrino factory of this kind, by comparing the oscillation probabilities of 
CP-conjugate channels z/j — >■ Vj and Vi — > Vj with [i ^ j). The most suitable channels 
for studying CP violation are Uf, — > z/^ and Ue —>■ i^fi, as well as their T conjugate 
partners z/^ v^. and z/^ v^- In these channels, unlike the channels, 
the CP violating part of the oscillation probability is not hidden by the CP conserving 
part 0, so that large asymmetries between CP-conjugate channels may arise, provided 
the leptonic CKM matrix allows for large violation of CP. Among these, the channels 
and z^e z/^ seem to be the most promising because it may be easier to 
disentangle negative from positive muons, in a large detector of high density, than to 
disentangle electrons from positrons. 

As a matter of example, we study the case of neutrino beams from stored muons 
with energy = 20 GeV, and an experimental arrangement where they travel over a 
distance of some 730 km, i.e. from CERN to the Gran Sasso laboratory, or, likewise, 
from FNAL to the Soudan mine. But we also give the scaling behaviour of the effects 
with either energy or distance. Within the range of squared mass differences and 
mixing angles, matter effects are important, see also However, unlike the case 
where neutrinos traverse the Earth from the antipode 0, they are easy to cope with 
because the density in the Earth's crust is essentially constant. 

The z/g oscillation probability can be measured as follows: Suppose the elec- 

tron neutrinos are produced by the decay of a number of positive muons in the 
straight section of the storage ring pointing to the detector. The which appear when 
there is oscillation of z/g into z/^, are detected by their charged current interaction in 
the detector. The number of observed muon neutrinos, n^^ is given by, 

n,,^ = 10^ Na j F,^ a,^ P{p, z/^) dE (1) 

where F^^ is the forward flux of electron neutrinos from a number A''^+ of positive 
muon decays, a^^^ is the charged current cross section per nucleon and P(z/e — > u^) is 
the oscillation probability for neutrinos traveling inside the Earth taking into account 
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matter effects. Nut is the size of the detector in kilotons. An analogous formula holds 
in the case of anti-neutrinos. 

Adopting the sample design configuration for muon production, cooling, accelera- 
tion and storage described by Geer 0, the number of available muons of either sign is 
approximately 8 ■ 10^° per year, for muons stored at an energy of 20 GeV. Of these, one 
fourth decay in a straight section directed towards the neutrino detector with a 10 kT 
target some ~ 730 km downstream, yielding about 2 ■ 10^° neutrinos per year and an 
identical number of anti-neutrinos. We use these numbers in what follows, and refer 
the reader to for details of the design of neutrino beams. 

Let us begin by computing the number of produced muon neutrinos. Experimental 
cuts needed to eliminate background as well as detecting efficiencies will be included 
later on. The neutrino fluxes at a neutrino factory have simple analytical forms that 
follow from the well-known formulae for muon decay. Let x = Ey/ be the fractional 
neutrino energy. For unpolarized muons of either charge, and neglecting corrections of 
order mf^/E'^ the normalized fluxes of forward moving electron neutrinos are 

^.^,p^(x) = 12x2(l-x) (2) 

and, for each neutrino type, the flux in the forward direction due to A*"^ decaying muons 
is 



dx dQ 



9u{x) (3) 



The above expressions are valid for a detector placed in the forward direction whose 
transverse dimensions are much smaller than the beam's transverse size ~ (L m^j/2 E^) . 

We assume that the interaction cross sections due to charged current interactions 
scale linearly with the energy even in the very low energy part of the spectrum 

a^, = .67 ■ 10~^^cm2 E (GeV) (4) 
ap^ = .34 ■ 10-^^cm2 E (GeV) (5) 

Regarding matter effects, let us remind the reader of the fact that of all neutrino 
species only z/g and Ue have elastic scattering amplitudes on electrons due to charged 
current interaction. This, as is well known, induces effective "masses" fi = ±2Eiya, 
where the upper sign refers to the electron neutrino, the lower sign to the corresponding 
anti- neutrino, and where a = y^GpUe, Ue being the electron density. 

Matter effects are important provided the interaction term fi, 



is comparable to, or bigger than, the quantity A„2 = rn^ — rn?^ for some mass difference 
and neutrino energy. 
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CP related observables often involve the comparison between measurements in two 
charge-conjugate modes of the factory. One example of an asymmetry is 

tot ^ JPj^e v^,)F^a^JE - JPjue Uf,)Fp^apJE 

/ P(z/e ^ iy^)F,^a,JE + / P(z7, ^ u^)F,^apJE ^ ' 

or in other terms, 

tot _ n.jN^+ - npjN^- 

In vacuum this quantity a^p would be a pure CP odd observable. The voyage through 
our CP uneven planet, however, induces a nonzero asymmetry even if CP is conserved, 
since z/g and Ug are affected differently by the electrons in the Earth [|1^. Therefore, to 
obtain the genuine CP odd quantity of interest, the matter effects must be subtracted 
with sufficient precision. 

For this purpose, we compute the matter asymmetry in the absence of CP violation, 
or fake CP asymmetry, by 



acp[<J = 0) = -pTTT r-j 77; — ; — . — r-j ^ wJ 





z^m) I<5=o Fy^a^^dE 


- 


z^m) I<5=o Fp^ayJE 




^'m) 1 5=0 Fy^ayJE 




i^fj) 1^=0 Fy^Gp^dE 



where we take into account matter effects but set 5 = in the transition probabilities. 

The total asymmetry a^p that will be found in an experiment of the type described 
above, is a function of acp{5 = 0), eq. (||), and of the asymmetry in vacuum (taking 
due account of CP violation) 

.ac ^ IP^'^'^i^e - y,) P.. g.. dE - /P-(^, ^ y,) Fp, a,^ dE 

/ P--(z/e ^ v^) Fy^ a,^ dE + J P--<^{ve ^ v^) Pp^ dp, dE ^ ' 

where P'"°-'^{i'f. —>■ v^) and P''"^(z/e — ^ v^j^) are the oscillation probabilities in vacuum. 
Provided acp{5 = 0) is not too large. 

a:;,^^ ^ a'S'p - acp{S = 0) (11) 

is a good approximation. In any case, the error one makes in calculating acp by means 
of eq. ([Tl|) is smaller than the uncertainties on acp{S = 0) itself. In addition, the error 
can be estimated by calculating the T-odd asymmetry for each neutrino energy, 

aT[Ey,d}-— r ■ — r [iZ) 

where a nonzero value cannot be induced by matter effects. This also means that 
a cleaner quantity in testing T violation than is acp for CP violation. 

An important component of any study of muon appearance due to z/g os- 
cillations is the event selection strategy for the /x's produced from charged current 
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interactions of the i^^'s. For neutrino experiments using a muon storage ring, the de- 
tailed prescription for event selection can be formulated only after the detector design 
is specified. There are, however, some basic issues concerning the signal and the back- 
grounds which all experiments are likely to be concerned with. On general grounds 
the background to a wrong sign muon signal is associated with the numerous decay 
processes that can produce fake muons: pions maskerading as muons, muonic charged 
currents (here one would also have to miss the right sign muon) or electronic charged 
currents (here one has also in the decay of the latter right sign muons), to name only a 
few. Without referring to a specific detector and the corresponding simulation toil we 
trust the experimental proficiency by setting an overall detection efficiency of 30% and 
by making a cut E^, > b GeV to eliminate inefficiently observed low energy interactions. 



2 Who mixes, two, three, or four flavours? 

With growing evidence for non vanishing neutrino masses, experimental studies of 
neutrino oscillations, and their analysis in terms of three (or more) flavours, have 
become popular and will continue to be of central significance for lepton physics in the 
future. 

The easiest way to describe any individual case of oscillations is to use a scheme 
where only two neutrino flavours are allowed to mix. Indeed, much work was done 
on analyses of neutrino oscillations in terms of two flavours but, as was pointed out 
by us and by others, the results for the squared mass differences may be misleading 
when applied to the real lepton world which contains three flavours. We summarize 
the situation regarding the squared mass differences as follows. 

In models involving three oscillating flavours one often relies on squared mass differ- 
ences which are taken from analyses of individual oscillation experiments in the frame- 
work of a two-flavour scheme. For instance, if one assumes the solar, the atmospheric 
and the LSND oscillations to be governed by just a single oscillation "frequency", 
AA4^, then the characteristic frequencies of the three oscillations, i.e. 

AA^Lio. = lO-^V^ or 10-W 



10-^eV^ < AAl|,^^3pi^^^i^ < 10-eV^ 

lO-^eV^ < AA^f^gj^Q < lO^eV^ 

cannot be reconciled with just three neutrino mass eigenstates. Therefore, in order to 
simultaneously accommodate all three oscillations as observed, under the assumption 
stated above, we must introduce (at least) a fourth neutrino. Since we know from the 
width of the boson that only three neutrino species have normal weak interactions, 
this extra, fourth neutrino must be sterile. 
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As there is no other, direct evidence for the existence of one or more sterile neutrinos, 
one is lead to conclude that assuming all observed oscillation phenomena to involve 
but a single AJ^^ is erroneous. Suppose, instead, that there are only three neutrinos, 
with masses such that 

ml-ml = AM^ < - = Am^ (13) 

Then, as was shown in it is possible to explain the LSND result as an oscillation 
involving AM^, the flavour conversion of solar neutrinos as one involving Am^ and 
the atmospheric neutrino anomaly as a mixture of both frequencies. In contrast to 
these findings, an analysis of the atmospheric data assuming (erroneously) that only 
one AA4^ is involved would find a value intermediate between those corresponding to 
the LSND and solar effects, as observed. 

In calculating observable effects of CP violation in neutrino oscillations we assume a 
scenario with three flavours, where the two squared mass differences obey the inequality 



(|T^) and lie in the range 

10"^ eV^ < Am^ < 10"^ eV^ , AM^ ^ 0.3 eV^ . (14) 

If there are three Dirac neutrino types, then the flavour eigenstates are related to 
the mass eigenstates by a 3 x 3 unitary matrix 

/ C12C13 S12C13 si3e"*'^\ 

U = -S12C23 - Ci2S23Sl3e*'' C12C23 - Si2S23Sl3e*'' S23C13 (15) 

I iS iS } 

\ -512523 — Ci2C23Si3e — £12^236 — S12C23S13 C23C13 / 

where Cij = cos 6ij and Sij = sin 6ij. If the neutrinos are Majorana particles, there are 
two extra phases, but these do not affect oscillations p2| . 

The transition probability in vacuum for a neutrino changing from i>i (ui) to Uj (uj) 
is given by the sum and the difference of CP-even and CP-odd pieces, respectively, 

P{Ui Vj) = Pcpi^i Vj) + PGP{l^i i^j) (16) 
P{Vi Vj) = Pcp{Vi Vj) - PQp{Vi Vj), (17) 

where 

Pcp{i^i Vj) = - 4ReJf^ sin^ A12 - 4Re4^ sin^ A23 - 4ReJ^'i sin^ A31, 

(18) 

Pof'i.^i 'i^j) = -8(JijJsin Ai2sinA23sin A31, 
with J the Jarlskog invariant and 

Jkh = UikUljUjhUli 

A,j = AmlL/AE (19) 



k 
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3 Results for CP asymmetries with three flavours 



As stated above we henceforth assume a scenario with three flavours of neutrinos 



characterized by the squared mass differences (|T4D and the strong mixing found in |jT3 
solution I, where 

^12^35.5°, ^^23 - 27.3°, ^^13^13.1°. (20) 

Although a detailed comparison might need further, refined analysis, this range 
of squared mass differences and the set of mixing angles (^) describes all observed 
neutrino anomalies in an overall and satisfactory manner. Here we show that this 
same set of parameters predicts a CP asymmetry which may well be large enough to 
be detectable with neutrino beams from a muon storage ring as described in sect. 1. 

We organize the discussion of our results as follows: We first present our main result 
for the expected asymmetry. Next, the role of matter effects is illustrated by some 
examples, followed by a comparison of CP violating asymmetries with time reversal 
violating asymmetries. We then turn to a comparison with previous results and show 
why we find asymmetries which are sizeably larger than the ones estimated previously. 
Finally, some remarks on efficiencies and statistical uncertainties are added. 

Let n be the number of muons, n the number of antimuons detected in one year's 
time in a detector placed at 732 km from the collider. Assuming for a moment an 
efficiency of 100%, we find the asymmetry {n — n)/{n + n) shown in Fig. 1, as a 
function of the neutrino energy Ey. Part (a) of the figure refers to the lower limit 
Am^ = 10~^ eV^, part (b) refers to the upper limit Am^ = 10~^ eV^ of (0). The 
solid line corresponds to setting 5 = (no CP violation), the dashed line shows the 
full asymmetry, assuming 5 = 7r/2. As a matter of example Fig. 2 shows the absolute 
numbers n and rl, obtained in one year of running, under the same assumptions as 
before and for Am^ = 10~^ eV^. 

The role of matter effects: Clearly, the smaller matter effects the cleaner the mea- 
surement of the effects of genuine CP violation will be. As the asymmetry due to 
charged- current matter interaction grows faster than the CP asymmetry, as a function 
of the baseline, intermediate distances between collider and neutrino detector are pre- 
ferred over long distances. We illustrate this observation quantitatively by defining the 
asymmetry 

acp[Eu,o)-— ^ , p.- ^ (21) 

and by calculating the ratio acp^E^.Q) /[acp{Ey,Ti /2) — acp{Ey,Q)], as a function of 
Ey. Fig. 3 shows this quantity for a baseline of 732 km, part (a), and a baseline of 
7332 km, part (b), corresponding to the distance from FNAL to the Gran Sasso. In the 
case of the very long baseline, the CP asymmetry is completely swamped by matter 
effects. 
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CP- versus T- asymmetry: We also computed the T-odd asymmetry (|T^) for the ex- 
ample Am^ = 10~^ eV^ and the shorter baseline L = 732 km and compared it to the 



CP asymmetry (^) from which the matter effects were subtracted. We found them to 
agree within reasonable limits, thus corroborating the approximation (|ll]). 

Comparison with previous results: The authors of ref. |0 who use the following sample 
set of parameters 

Aml^ = 10"^ eV^ Amj^ = lO'^ eV^ sin^ 6^ = -5, ^23 = 45°, ^13 = 13° , (22) 

find CP violating effects which are markedly smaller than the ones we showed above. 
The explanation for this difference is simple: The most noticeable difference between 
the set ( P^ ) and ours is the value of Am"^^. Assuming the values ( P^ both the 
(12)- and the (13)-channels are strongly affected by matter effects, the effective sines 
{sin 9 12) matter and {sin 9 is) matter ^s defined in quickly tend to 1 as the parameter /x, 
eq. (|), increases. Consequently, the simultaneous interplay of all three fiavours and, 
hence, the visibility of CP violation decrease. In contrast to this situation our values of 
squared mass differences are such that only {sin 9 12) matter is affected while (sin 9i3)matter 
and, of course, {sin 923) matter remain unaffected. It is convenient to define an effective 
mixing matrix Vik which is obtained from (^) by replacing the sines and cosines by 
the matter affected sines and cosines, {sin 9 matters 

etc. 

To illustrate the comparison Fig. 4(a) shows the pertinent, effective matrix elements 
for our set of parameters, eqs. (|1^) and (|20D, as a function of /x, while Fig. 4(b) shows 



the same matrix elements for the set (p2D . In the latter case, both Vii and V12 tend to 
zero with /i increasing to its value in the Earth's crust. 

Efficiencies and statistical uncertainties: The discussion of statistical uncertainties 
is straightforward. In order to exclude the possibility that a measurement of a non 
vanishing genuine CP violation is due to a statistical fluctuation, the measured value 
must be larger than n ■ S{acp)stat, where S{acp)i,t3^t is the la statistical error on acp in 
the absence of CP violation, and n is the number of standard deviations we require in 
order to be happy with our result. Since in absence of CP violation the expectations 
of n,j^/N^+ and n^/N^- are equal, we get 

where (rij^^) {{np^)) is the expected number of z/^ (t'^) interactions seen in the detector. 

Regarding the contribution of the background to the statistical error, and according 
to the current estimates, the main source of background will be due to charm production 



in the charged current neutrino interactions in the detector 



—>■ i^fi ^ CC interaction —>■ fi 
c — > c decay 



found 
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Clearly, this background only affects the signal corresponding to the Ve ^^^i oscilla- 
tions (because we expect there /i+ appearance), and therefore, this background "noise" 
should be subtracted appropriately in the counting of rip^ in Eq.(||). Such a subtraction 
introduces a further source of statistical error. Using the estimate that 



rip, Lck^ 10 'rip^ |p=i (24) 

where np^ \p=i is the number of antimuon neutrino interactions that would be seen if 
all the initial antielectron neutrinos oscillated into antimuon neutrinos, we find 

as the complete expression for the statistical error. It is important to notice that even 
with only one year of data taking and a modest 30% detecting efficiency, the statistical 
error will be small enough to rule out the possibility of attributing to a statistical 
fluctuation a measurement of a non vanishing CP violation. 



4 Conclusions 

In summary, the two rather different "frequencies" (0), together with the strong mix- 
ing of all three flavours that describe the solar neutrino deficit, the atmospheric oscilla- 
tions, and the LSND anomaly, lead to relatively large CP and T violating asymmetries 
in neutrino oscillations. With the set of parameters (|l^ and ([20| ) the full interference 
of all three flavours is well developed and is only moderately damped by matter ef- 
fects. We have also tried the other solutions to the neutrino anomalies we had found 
in |13[ but find no more than 20% changes in the asymmetries. Among these the CP 
asymmetry seems large enough to be measurable with neutrino beams from a 20 GeV 
muon storage ring and with a detector at some 730 km from the source, corresponding 
to the distance of the Gran Sasso laboratory from CERN or, likewise, of the Soudan 
mine from FNAL. 
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Figure 1: (n — n)/{n + n) as a function of the neutrino energy for Am^ = 10"'^ eV^ 
(a) and Am^ = 10^'' eV^ (b). The sohd hne correspond to 5 = while the dashed hne 
correspond to 5 = 7r/2 
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Figure 2: Absolute number of muons and antimuons detected in one year's time 
in a 732 km baseline for Am^ = 10~^ eV^ and assuming 100% detecting efficiency, 
dashed and solid lines as before. 





Figure 3: "Noise over signal" ratio as a function of the neutrino energy for Am^ = 
10~^ eV^ and a detector placed at 732 km (a) and 7332 km (b) from the collider. 
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Figure 4: Effective matrix elements as a function of tlie effective mass for tlie set of 
parameters used tfirougli tliis work (a) and tliose of ref. @] 
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